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TECHNICAL FIELD 



5 The present invention relates to the automatics and computation techniques 

and, more particularly, to methods of the prehminaiy image processing for 
sharpening and contrast enhancement. 

BACKGROUND ART 

10 

One of the methods known in the prior art is described in the Russian 
Patent No.20 15561, pubhshed on 16.04.91, Int. Class G06K 9/40. According to 
this invention, the image correction is made basing on the analysis of the original 
image at the processing pixel and the local average value over some neighborhood 
15 of this pixel. 

The method of that patent smoothes the original image, thus producing the 
smoothed image containing the low frequency components presenting the image 
background. Then the smoothed image is subtracted from the original one 
producing the second image containing hi^ frequency components without 
20 background, said second image is then emphasized and added to the smoothed 
image. 
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The disadvantage of this method is that it emphasizes not only the useful 
signal but also the noise containing in high frequency image components, thus 

degrading the quality of the enhanced image. 

The method according to U.S. Patent No. 5,038,388, pubhshed on 
5 06.08.91, Int. Class G06K 9/40 smoothes the original image and subtracts the 
smoothed image from the original one thus producing the second image containing 
the high frequency image components only. Said second image is then adaptively 
5 emphasized so that the scaling factor is then higher tiian larger is a difference 

m between the processing pixels and its nei^borhood. The output image is 

J; 10 produced by adding the adaptively emphasized second image to the original 
%j image, thus sharpening the image without the noise amplification. 

y The solution described in this patent does not provide any noise 

13 suppression as tihis method can only emphasize the high frequency image 

^3 components that may contain a noise. 

1 5 Furthermore, the disadvantage of this method known in the art is that it fails 

to improve considerably the sharpness of weak edges as such enhaucement 
requires to emphasize the high frequency image components in the regions where 
a difference between the processing pixel and its neighborhood is comparable to 
the noise level. Therefore the edge enhancement in such regions causes the noise 

20 emphasis. 

According to the image processing method disclosed in U.S. Patent No. 
5,381,490, pubhshed on 10.01.95, Int. Class G06K 9/40 the largest difference A 



AMENDED SHEET 



-21-08-2000 RU 009900272 



PCT/RU99/00272 
Replacement sheet 

between the processing pixel and its nearest neighbors is calculated. Depending 
on the magnitude of this difference, one of the three processing modes is selected: 

- the edge enhancement by means of emphasis of the high frequency image 
components if A > Tj, here Tj presents the first pre-defined threshold value; 

5 - reproduction of the original image, if T2 < A < Ti, where T2 stands for the 
second pre-defined threshold value. 

- the image smoothing to si5>press a noise, if A < T2. 

The first disadvantage of this method is that it may emphasize the noise 
selectively if the difference A varies around any of the threshold values for some 
10 neighboring pixels thus transforming small differences between neighboring pixels 
into larger ones by involving different processing modes for these nei^boring 
pixels. 

Furthermore, this method fails to provide quahty enhancement of images with 
different noise magnitudes without tuning as the threshold values Ti and T2 are 
15 not selected adaptively. 

Another approach to the noise suppression in images is described in U.S. 
Patent 5,563,963, published on 08.10.96, Int. Class G06K 9/40. The method of 
this patent operates by selecting a plurality of groups of neighbors for each pixel 
of the original image, each group being a square containing N x N pixels, here N 
varies for example from 2 to 21. The processing pixel may be located at any 
position inside this square group of pixels. The least square best fit plane (the 
planar approximation) is used then to ^proximate pixel values in each of selected 
pixel groups and the new value for the processing pixel and the goodness of fit are 
computed basing on said approximation for each group. 



20 
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The target pixel of an enhanced image is produced by the weighted summing 
of all the new pixel values, the wei^t of a group is then higher than higher is the 
goodness of fit for this group. 

The disadvantage of this method is that it fails to enhance edges as this method 
5 provides the noise smoothing only. Furthermore, this method requires huge 
computation efforts to bufld least square approximations by hundreds of groups 
for each of hundreds of thousands of pixels. 

A technique for making low contrast parts of a digital image more 
perceptible is disclosed in the European Patent EP 0 756 247 Al, pubhshed on 
10 29.01.97, Int. CI. G06T 5/20. This method decomposes the original image into 
pluraUty of images in such a way that these images belong to different isotropic 
frequency bands. Thereat the image presenting the low frequency components of 
the original image has lower resolution (smaUer size) than the image presenting 
the high frequency components. The emphasis processing is carried out on these 
15 images. The degree of emphasis is selected to be lower for pixels having values 
larger than a predeteraiined direshold value than for pixels having smaU values. 
The output image is assembled from the processed images. 

The disadvantage of this method is that it emphasizes noise as weU as the 
useful signal. This method fails to provide the edge detection and enhancement as 
20 the isotropic frequency channels are used. 

The method disclosed in U.S. Patent 5,739,922, pubhshed on 14.04.98, Lit. 
Classes G06K 9/40, H04N 1/40 operates by sphtting of the original color image 
into three isotropic frequency channels: the low frequency image components 
(LF), the medium frequency components (MF) and the high frequency 
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components (HF). The adaptive emphasis of the HF components and adaptive 
suppression of the MF components is then carried out, thereat the multiphers for 
HF and MF image components are then higher than higher is a correlation 
between at least two of three basic image colors. The enhanced image is obtained 
5 by s ummin g LF image components with adaptively suppressed MF components 
and adaptively emphasized HF image components. 

The image processing method and apparatus described in said patent may have 
the hmited appHcation as they are suitable for color images only as the correlation 
between color components is used for carrying out the image processing. 
10 Furthermore, the noise suppression according to this invention is significantly 
hmited, as the HF image components, that also contain a noise, may be 
emphasized only and the noise suppression in MF image components is hmited 
because no directional sphtting of the original image is used. 

The edge detection and enhancement can not be obtained by this method as the 
15 isotropic frequency chaimels are used. 

All these disadvantages degrade the quahty of enhanced images. 

The most relevant image processtag method is described in U.S. Patent No. 
5,351,305, pubhshed on 27.09.94, Int. Class G06K 9/40. According to this patent, 
a plurahty of directionafly filtered images is obtained fi-om the original image by 
20 applying directional filters in the fi-equency domain. The enhanced image is then 
formed by selecting each target pixel either from the directionaUy filtered image, if 
the contrast edge is detected nearby the processing pixel or from the original 
image otherwise. Thereat the contrast edge is detected nearby the processing pixel 
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by generating the standard deviation image and by producing the eigenvector 
description of this image. The eigenvector length is compared to the pre- 
determined threshold value to detect an edge. 

The target pixel is equal to the corresponding pixel of the original image, if the 
5 edge was not detected nearby. Otherwise, the target pixel is selected from an 
image filtered with the most nearly corresponding direction of filtering. 

While detecting edges, the eigenvector length may vary around the threshold 
value for several adjacent pixels. Thereby die neighboring pixels of the enhanced 
image are selected from different images (the original image and directionaUy 
10 filtered image) thus causing the selective noise emphasis. This emphasis degrades 
the enhanced image quahty. 

Furthermore, original images may differ in their noise levels thus requiring 
different threshold values. This method does not include the adaptive selection of 
the threshold value and therefore may not provide high quality processing of 
1 5 images with different noise levels. 

Provided that the edge is detected nearby, the selection of pixels of the 
enhanced image is made from one of plurahty of directionaUy filtered images thus 
causing the complete suppression of aH image structures that differ by their 
direction firom the detected edge, notwithstanding fliat those structiires can be 
20 clearly seen in the original image. 
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DISCLOSURE OF INVENTION 

The object of the claimed invention is to provide a method for enhancing 
the image sharpness and contrast combined with simultaneous noise suppression. 
5 The said objective is achieved in the image processing method, comprising 

the presentation of the original image as a matrix of the discreet picture elements 
(pixels), sphtting of the original image into n frequency channels, each channel 
presented by an image matrix of the same size as the original image, the edge 
5 detection and assembhng of the output image from said n frequency channels 

I 10 taking the detected edges into account, by sphtting of the original image into the 
; j: low frequency channel and n-1 high frequency channels and edge detection by 

computation in each of n-1 selected high frequency channels of the correlation 
■ between processing pixel and its neighboring pixels, followed by the comparison 
of said correlation vahie with that for the corresponding pixels in other high 
15 frequency channels and with the threshold vahie for this channel. Basing on the 
results of comparison, the weighting coefficients are formed for each pixel of each 
of n-1 high frequency channels, and the assembly of the output image is made by 
s ummin g each pixel from the low frequency channel with all products of the 
corresponding (by their location in the image) pixels of n-1 high frequency 
20 channels by their wei^ting coefiBcients. 

The said objective is also achieved by selection of m of n-1 high frequency 
channels (2<m<=n-l) in such a way that they differ one from another in die 
direction of principal passing only. Therewithal, the weighting coefficients for any 
of pixels of any of m high frequency channels is defined basing on the comparison 
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of its correlation value to the threshold value and to the correlation values of the 
corresponding (by their location in the image) pixels of other m-1 high frequency 
channels. 

The said objective is also achieved by computation of said correlation as a 
5 product of the processing pixel value and the v/eighted sum of its neighbors, said 
weights being anisotropic thereat the direction of this anisotropy corresponds to 
the direction of principal passing of the processing channel. 

1 Furthermore, the threshold value for any of n-1 high frequency channels is 
.| determined by analysis of distribution of values, or absolute values, of pixels in 
■y 10 this channel. The threshold value may be also determined by analysis of 

2 distribution of values, or absolute values, of pixels of the original image. 

□ The claimed as the invention image processing method is characterized by the 

following features that distinguish it from the most relevant method known in the 
1:3 prior art: 

15 1. while splitting the original image into frequency channels, the low frequency 
channel and n-1 high frequency channels are selected; 
2. the value of correlation between the processing pixel and its nei^boring pixels 
is used to detect edges in the image. More specifically, the value of said 
correlation is compared to the correlation values for the corresponding (by 

20 their location in the image) pixels in other n-1 high frequency channels and to 
the threshold value for this channel. The weighting coefBcients for all pixels of 
n-1 high frequency channels are generated as a result of this comparison. 
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3. the oix^ut image is assembled by means of summing of each pixel jfrom the 
low frequency channel and all the corresponding (by their location) pixels from 
n-1 high frequency channels multiplied by their weighting coefficients. 
Concerning the first feature it should be noted that the extraction of the low 
frequency channel that is not a subject of any processing provides the distortion 
free reproduction of large-scale objects of the original image in die enhanced 
image, as these objects are passed throu^ the low frequency channel without any 
processing. Therewithal the energy of noise passing through the low frequency 
channel without siippression is low as most of the noise energy is concentrated at 
high frequencies. 

Furthennore, the extraction of the low frequency channel helps reduce the 
computation efforts needed to find correlation values for n-1 high frequency 
channels, as the subtraction of the local average value is required to compute 
correlation. This subtraction is made while extracting (subtracting from the 
original image) the low frequency channel. 

The spHtting of the original image into several high frequency channels helps 
improve, compared to selection of one or two channels, the noise suppression, as 
a noise associated with pixels of those frequency channels where edges were not 
found nearby, is prevented from contributing to the enhanced image. 

Concerning the second feature it should be noted that edge detection by means 
of the correlation value between the processing pixel and its neighbors allows to 
find weak edges on a noisy background as the correlation value is close to zero 
for the typical noise, thereat the correlation values for adjacent pixels forming the 
edge are positive and of similar magnitude. This difference in correlation behavior 
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helps achieve the high noise suppression combined with edge emphasis thus 
increasing the image sharpness and contrast. 

The edge detection basing on the correlation between a processing pixel and 
its neighbors provides a quantitative measure of tbe edge intensity for each image 
pixel. This allows performing the selective emphasis of edges having certain 
intensity, for example, the weak edges. Furthermore, use of the correlation 
between the pixel and its neighbors makes the method of the present invention 
applicable to various types of images, including, for example, color and gray scale 
images, as well as three-dimensional hnages. 

The third speciJQc feature, namely assembling of the enhanced image by 
wei^ted summing of images from all frequency channels, helps to remove 
completely the effect of the selective noise emphasis. This effect is caused by 
selection of neighboring pixels from different frequency channels. 

Furthennore, the determination of the threshold value by means of analysis of 
statistical distribution of pixels provides a high quahty processing of images with 
materially different noise magnitudes without changing of parameters. 

Fmthermore, use of anisotropic frequency channels and anisotropic wei^ts 
makes the image processing method highly sensitive to weak edges. 

BRIEF DESCRIPTION OF DRAWINGS i 

i 

The claimed invention is illustrated by the drawings of the apparatus that 
implements the claimed method. 
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Fig.l is a block diagram of the apparatus. 

The preferred embodiment of the sub-units of said apparatus are shown in 
more details at Fig. 2-5: 

Kg.2 is a block diagram of the frequency channel splitting unit (sphtting unit). 

Hg3^s a block diagram of the unit for computation of correlations (correlation 
unit). 

Fig.4 is a block diagram of the unit for forming of the weighting coefficients 
(wei^ting imit). 

Fig£is a block diagram of the unit to assemble the output image (assembhng 
unit). 

Fig.6 shows an example of the pre-defined channel selection matrixes. 

^jTiUustrates the operation of the frequency channel splitting unit. 

Fij^s a graph showing an example of the dependence of the weighting 
coefficient on the correlation value. 

Referring to Fig.l, the apparatus contains the image source 1, the output of 
said image source connected to the input of the splitting unit 2. The low frequency 
output 7 of said unit 2 is connected to the input 9 of the assembhng unit 5, all 
other outputs of said spUtting unit are connected to the corresponding inputs of the 
correlation unit 3, said other outputs are also connected to the inputs lOi - IO4 of 
the assembhng unit 5. The outputs of said correlation unit 3 are connected to the 
corresponding inputs of the weighting unit 4, the outputs of said weighting unit are 
the inputs 81 - 84 of the assembhng unit 5. Thereat the output of said assembling 
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unit 5 is connected to the input of the memory unit 6, the output of said memory 
unit is an output of the apparatus. 

Fig.2 shows more in details the preferred embodiment of the splitting unit 2. 
It includes the direct Fourier processor 11, the input of said processor is 
5 connected to the output of image source 1, the output of said direct Fourier 
processor 11 is coimected to the first inputs of matrix multipKcators 12o - 124 The 
second inputs of said multiplicators are connected to the corresponding memory 
i;3 units, said memory units hold the pre-defined channel selection matrixes 

g iUustrated at Fig.6. Any of the matrix multiplicators 12o - I24 performs the 

10 element-by-element multiplication of matrixes supplied to its two inputs. The 
outputs of said matrix multiplicators 12o - I24 are connected to the inputs of the 
inverse Fourier processors 14o - I44. The output of the inverse Fourier processor 
Q 14o is connected to the input 9 of the assembling unit 5, and outputs of the inverse 

m Fourier processors 14i - 1 44 are connected to the inputs lOi - IO4 of the 

p 15 assembling unit as weU as to the corresponding inputs ofthe correlation unit 3. 

Fig.3 shows more in details the correlation unit 3. The memory unit 15 
holds the image of the processiug frequency channel. The input of said memory 
unit is the input of the correlation unit. It is connected also to the input of the 
noise level measuring unit 20, the output of said unit is connected to the first input 
20 of divider 19. The second input of said divider is connected to the output of 
multiplicator 18. The first input of said multipHcator is connected to the first 
output ofthe memory unit 15. Other outputs of said memory unit are connected to 
the weighting summator 17, the output of said weighting summator is the second 
input of multipHcator 18. The address input ofthe memory unit 15 is connected to 
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the address generator 16. The output of divider 19 is the output of the correlation 
unit. 

The noise level measuring unit 20 may be implemented according to the 
U.S. Patent No. 5,657,401, published on 12.08.97, Int. Qass G06K 9/40. 

5 AH the memory units are of the random access memory type and they are 

weU known in the art. 

The weighting summator 17 may be implemented as eight scalar 

O 

a multiphcators (the number of said muItipHcators is equal to the number of 

i:| adjacent pixels). Any of said scalar multipHcators has two equal inputs and one 

=p 10 output. The outputs of all said multiplicators are connected to the inputs of 
summator, the output of said summator is the output of unit 17. The first inputs of 
Q said multiphcators are the inputs of the unit 17 and the pre-defined weighting 

Q coefficients are suppKed to the second inputs of said multiplicators. 

y Fig-4 shows the weighting unit 4. The four inputs of the weighting unit 4 

15 are the inputs of four rounding imits 23 1 - 234, the outputs of said rounding units 
are connected to the inputs of the address assembling unit 24. The output of said 
unit 24 is connected to inputs of four look-up tables 25i - 254. The look-up table is 
a memory unit that stores the values of the weighting coefficient for any set of the 
four input correlation values. The outputs of the look-up tables 25i - 254 are 
20 connected to the inputs of the memory units 26i - 264, said memory units 
accumulate values of weighting coefficients. TTie address inputs of the memory 
imits 26i - 264 are connected to the address generator 27 and the outputs of said 
memory units are connected to the inputs of summators 28 1 - 284 for averaging of 
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weighting coeflScients. The outputs of said sunimators are the outputs of the 
weighting unit 4. 

Fig.5 shows the assembling unit 5. It consists of four multiphcators 29i - 
294 and summator 30. The first inputs 8i - 84 of said multiphcators are connected 
to the outputs of the weighting unit 4 and the second inputs lOi - IO4 of said 
multiphcators are connected to the outputs of the splitting unit 2. The outputs of 
multiplicators 29i - 294 are connected to the corresponding inputs of the 
summator 30 and the input 9 of said summator is connected to the output 7 of the 
sphtting unit 2. The output of summator 30 is connected to the input of the 
mraiory unit 6 that accumulates an enhanced image. 



BEST MODE FOR CARRYING OUT 
THE BSrVENTION 

The apparatus implements the claimed method as it is described hereinafter. 
Referring to Fig.l, the input image is generated by the image source 1. The MRI 
unit may be used, for example, as an image source 1 . Said MRI unit produces the 
image of the oross-section of an object, this image being a matrix containing 
discreet picture elements (pixels). Said image is carried to the input of the splitting 
unit 2. The operation of the splitting unit 2 is described with a reference to Fig.2 
and Fig.7. The input image is transformed to the frequency presentation by the 
direct Fourier processor 11. Said frequency presentation contains the conq)lete 
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information about the original image and is represented by the matrix of the same 
size as an input image. This matrix is passed to identical matrix multiplicators 

12o - 124, these matrix multiphcators perfbnn the element-by-element 
multiplication of the frequency presentation of the original image by pre-defined 
5 channel selection matrixes. Said channel selection matrixes are stored in memory 
units 13o - 134. Each channel selection matrix contains the multipliers for all 
spatial frequencies of the frequency image presentation. Fig.6 shows the examples 
of the channel selection matrixes. 

More specifically, as an image is presented by a 2D matrix, its frequency 
10 presentation is also a 2D matrix. Fig.6a shows schematically a frequency 
presentation matrix. The horizontal and vertical spatial frequencies vary along 
axes kx and ky, respectively. 

The zero spatial frequency corresponds to the constant image density. It is 
located at the crossing point (31) of axes kx and ky. 

15 Points 32 and 33 represent the largest spatial frequency in horizontal 

direction. The examples of images contributing to these points are shown in 
drawings 34 and 35. 

Similarly, the maximal spatial frequency in the vertical direction is located 
at points 36 and 37; the example of image contributing to these points is 
20 illustrated by drawing 38. 

The maximal spatial frequencies are located at points 39-42. The example 
of an image contributing to these maximal spatial frequencies is shown in drawing 
43. 
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Medium spatial frequency in horizontal direction is located at point 45. The 
example of image contributing to this point is shown in drawing 44. 

The location of the spatial frequencies in drawings Fig.6(b-f) corresponds 
to the scheme depicted in Fig.6a. 

5 Fig.6b shows schematically the pre-defined selection matrix for the low 

frequency channel, said matrix being stored in the memory unit 13o. 

The dark area 46 is jSIled by the unit values of the matrix elements. This 
area corresponds to spatial frequencies that pass through the low frequency 
channel. The white region is filled by the zero values of the matrix elements, 
10 therefore the frequencies of the white region do not pass through the low 
frequency channel. 

Fig.6(c-f) show schematically the selection matrixes for four high frequency 
channels, thereat the same notations as in Fig. 6b are used. 

It should be noted that the sum of all channel selection matrixes Fig.6(b-f) 
15 is the matrix with all elements equal to 1. Therefore all the information from the 
original image passes through at least one channel. 

Referring to Fig.2, each of the matrix multipUcators 12o - I24 forms on its 
output the matrix of the corresponding frequency channel in the frequency 
presentation. The inverse Fourier processors 14o - 144 transform said matrixes to 
20 the coordinate presentation. 

The direct Fourier processor 1 1 and inverse Fourier processors 14o - 144 
may be implemented basing on the Fast Fourier Transform algorithm as described. 
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for example, in: Cooley, J.M., Lewis, P.A.W. and Welch, P.D. (1969) The finite 
Fourier transform. I.E.E.E. Trans. Audio Electroacoustics AU-17, 2, 77-86. 

Fig.7 further illustrates operation of the splitting unit. Fig.7a shows the 
example of an input image, Fig.7(b-f) show the images formed on outputs of the 

5 inverse Fourier processors 14o - 144, respectively, as a result of processing of the 
image shown in Fig.7a. The image of a low frequency channel 7b is carried from 
the output of the Fourier processor 14o, being the output 7 of the splitting unit 2, 
to the input 9 of the assembly unit. The images of four high frequency channels 
are carried from the outputs of the Fourier {Mrocessors 14i - I44, being another 

0 outputs of the sphtting unit 2, to the coixesponding inputs of the correlation unit 3 
and to the inputs lOj - IO4 of the assembly unit 5. 

The further processing of these images will be described on example of one 
high frequency channel as said processing is identical in aU high frequency 



15 Referring to Fig.3, the memory unit 15 stores the partial image of the 

processing channel. To compute the unnormahzed correlation value, the center 
pixel value 21 and values of its neighboring pixels 22 are sequentially selected 
from the memory 15. Said values of neighboring pixels pass to the input of 
weighting summator 17, said summator implements the foUowing operation on 

20 pixel values: 



channels. 



N 



r = 



i = l 
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where N is a number of pixels in neighborhood 22 of the central pixel (preferably 
N=8), Vi are the pre-defined weighting coefiScients (preferably Vi=l/8) and Xi are 
the pixel values fi-om neighboriiood 22. The multipUcator 18 forms a product of 
the wei^ted sum of neighboring pixels and the central pixel value. This product is 
the unnormalized correlation value for the processing pixel. It is compared to the 
threshold value by dividing by said threshold value (output of the noise level 
measuring means 20) in divider 19. The result of said division is compared to 1.0 
in the weighting unit 4. The processing described herein is repeated for all pixels 
of the partial image of the processing frequency channel. 

The image of this fi-equency channel passes also to the noise level 
measuring means 20. The noise level from the output of means 20 is used as a 
threshold value to normalize correlation values by divider 19. As a result, the 
matrix containing the correlation values for all pixels of the processing firequency 
channel is formed on output of the correlation unit 3, said correlation values being 
normalized by the threshold value for the processing frequency channel. 

The correlation values formed by the correlation unit 3 are carried to the 
weighting unit 4. Referencing to Fig.4, said correlation values for four frequency 
channels pass to inputs of the rounding means 23 1 - 234. Said rounding means 
decrease the data precision to 4 or 5 bits. 

The four rounded values from outputs of means 23 1 - 234, each containing 
4 or 5 bits, are assembled into one 16 or 20 bit word by the address assembling 
unit 24. The formed address is used as an input value for four look up tables 25i - 
254. Each of said look up tables is a memory unit that stores values of the 
weighting coefScients for any combination of four correlation values in four 
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frequency channels, thereat such combination defines the address formed by the 
unit 24 in a unique way. 

Fig. 8 shows the preferred dependence of the weighting coefficient Wi in any of 
the frequency channels on the correlation value Q in this channel and correlation 
5 values in other three channels, where A stands for the threshold value for this 
frequency channel. The weight Wi depends on the correlation value Q and 
maximal correlation value L in other three frequency channels. This dependence is 
illustrated at Fig.8 by the curves: 

H curve A for Q > 0.7 L, 

m 10 curve B for Ci = 0.5 L, 
H curve C for Ci = 0.3 L, 

ii curve D for Q = 0. 1 L and 

W curve E for Ci = 0.01 L. 

[U The memory units 26i - 264 accumulate values of the weighting coefficients 

u 15 generated by look up tables 25 1 - 254. The address generator 27 and summators 
28i - 284 smooth said weighting coefficients in each frequency channel. Said 
smoothing is obtained by summing in the summator (for example, 28i) the center 
value of the coefficient and its neighboring values being sequentially selected from 
the memory unit (for example, 26 1) by the address generator 27. The smoothed 
20 values of weighting coefficients formed on outputs of summators (for example, 
28i) pass to the outputs of the wei^ting unit 4. 

The operation of the assembling unit is described with a reference to Fig.5. 
The values of weighting coefficients for four frequency channels pass from the 
outputs of the unit 4 to the inputs 81 - 84 of multiphcators 29i - 294 and the pixel 
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values of the corresponding frequency channels are carried from outputs of the 
spHtting unit 2 to another inputs 10] - IO4 of said multipHcators. The products of 
said pixel values by the corresponding weighting coefiBcients generated by 
multiplicators pass to the inputs of the summator 30. Thereto the corresponding 
5 pixel value of the low frequency channel passes to the input 9 of said summator. 
Thereby the summator 30 calculates the value of the target pixel of the output 
image, said pixels are accumulated in the memory unit 6. 

The embodiment described herein ill^Ist^ates the method as applied to 2D 
scalar images. It is understood however that the claimed method may be appUed 



X; 10 similarly to 3D images, hi this case, in the apparatus used to implement the 
claimed method: the number of frequency channel increases, the 3D Fourier 



(used, for example, to compute a correlation value) is 26 instead of 8. 

The claimed method may be apphed also to processing of vector images, 
15 particularly the color images. Thereat, the 3 components of a vector presenting a 
pixel value may correspond for example to the intensity of the 3 basic colors for 
this pixel. In this case, the scalar operations on pixel values, like Fourier transform 
and summing, are replaced by the corresponding vector operations as it is known 
in the art and the correlation is computed as a scalar product of the center pixel 
20 value and the weighted vector sum of its neighbors, thereat the vector summator 
contains as many scalar summators as the number of vector components. 



processors are used iustead of 2D ones and the number of neighboring pixels 
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INDUSTRIAL APPLICABILrrY 



The image processing method according to the invention has the following 
advantages. 



neighboring pixels helps detect weak edges against a noisy background. This 
feature provides high noise suppression in conjunction with the emphasis of weak 
edges, thus significantly improving the image quality. 

Furthermore, use of correlation between the processing pixel and its 
10 neighbors makes the method appHcable to a wide variety of image types, 
including color images, gray scale images and 3D images. 

Second, the determination of the threshold value by analysis of distribution 
of pixel values provides a high quality processing of images with the material 
difference in their noise levels without changing of parameters. 

15 Third, the distortion-fi-ee reproduction of the large-scale image structures is 

achieved due to the separation of the low frequency channel. 



5 



First, use of the correlation between the processing pixel and its 
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